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Abstract 

In illiquid markets, option traders may have an incentive to increase their portfolio value 
by using their impact on the dynamics of the underlying. We provide a mathematical frame- 
work within which to value derivatives under market impact in a multi-player framework by 
introducing strategic interactions into the Almgren & Chriss (2001) model. Specifically, we 
consider a financial market model with several strategically interacting players that hold Eu- 
ropean contingent claims and whose trading decisions have an impact on the price evolution of 
the underlying. We establish existence and uniqueness of equilibrium results and show that the 
equilibrium dynamics can be characterized in terms of a coupled system of possibly non-linear 
PDEs. For the linear cost function used in Almgren & Chriss (2001), we obtain (semi) closed 
form solutions for risk neutral or CARA investors. Finally, we indicate how spread crossing 
costs discourage market manipulation. 
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1 Introduction 



Standard financial market models assume that asset prices follow an exogenous stochastic process 
and that all transactions can be settled at the prevailing price without any impact on market 
dynamics. The assumption that all trades can be carried out at exogenously given prices is 
appropriate for small investors that trade only a negligible proportion of the overall daily trading 
volume; it is not appropriate for institutional investors trading large blocks of shares over a short 
time span. The trading decisions of institutional investors are likely to move stock prices in an 
unfavorable direction and often induce significant trading costs. 

It is now widely recognized that (the lack of) liquidity is a major source of financial risk 
and there has been an increasing interest in mathematical models of illiquid financial markets. 
Much of the literature on illiquidity focusses on either optimal hedging and portfolio liquidation 
strategies for a single large investor under market impact (Cetin, Jarrow & Protter (2004), Alfonsi, 
Fruth Sz Schied (2007), Rogers Sz Singh (2007)), predatory trading (Carlin, Lobo & Viswanathan 
(2007), Schoeneborn & Schied (2008)) and the role of derivative securities including the problem 
of market manipulation using options (Jarrow (1992), Kumar & Seppi (1992)). In an illiquid 
market, derivative traders have an incentive to utilize their impact on the price dynamics of the 
underlying in order move the option value in a favorable direction^} It has been shown by Jarrow 
(1994), for instance, that by introducing derivatives into an otherwise complete and arbitrage-free 
market, manipulation strategies with a risk free gain may appear, such as market corners and front 
runs. Schoenbucher & Wilmott (2000) discuss an illiquid market model where a large trader can 
influence the stock price with vanishing costs and risk. They argue that the risk of manipulation 
on the part of the large trader makes the small traders unwilling to trade derivatives any more. In 
particular, they predict that the option market breaks down. Our analysis indicates that markets 
do not necessarily break down when stock price manipulation is costly as it is in our model. 

While the aforementioned papers differ significantly in their degree of complexity, they all focus 
on a single player framework. When multiple players are considered the analysis is typically con- 
fined to some form of "stealth trading" as in Carlin, Lobo & Viswanathan (2007) and Schoeneborn 
& Schied (2008), where liquidity providers try to benefit from the liquidity demand that comes 
from some "large" investor but no strategic interaction between liquidity suppliers and consumers 
is considered. In fact, so far only little work has been devoted to models with strategically in- 
teracting market participants. Vanden (2005) considers a pricing game in continuous time where 
the option issuer controls the volatility of the underlying but does not incur liquidity or spread 
crossing costs. He derives a Nash equilibrium in the two player, risk neutral case and shows that 
"seemingly harmless derivatives, such as ordinary bull spreads, offer incentives for manipulation 
that are identical to those offered by digital options" (p. 1892, 1. 36). Gallmeyer & Seppi (2000) 

1 Gallmeyer & Seppi (2000) provide some evidence that in illiquid markets option traders are in fact able to 
increase a derivative's value by moving the price of the underlying. 
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consider a binomial model with three periods and finitely many risk neutral agents holding call 
options on an illiquid underlying. Assuming a linear permanent price impact and linear transac- 
tion costs, and assuming that all agents are initially endowed with the same derivative they prove 
the existence of Nash equilibrium trading strategy and indicate how market manipulation can be 
reduced. 

We provide a general mathematical framework within which to value derivative securities in 
illiquid markets under strategic interactions thereby extending the work of Gallmeyer & Seppi 

(2000) in several directions. Specifically, we consider a pricing game between a finite number of 
large investors ("players") holding European claims written on an illiquid stock. Their goal is 
to maximize expected utility at maturity from trading the stock where their portfolio value at 
maturity depends on the trading strategies of all the other players, due to their impact on the 
dynamics of the underlying. Following Almgren & Chriss (2001) we assume that the players have 
a permanent impact on stock prices and that all trades are settled at the prevailing market price 
plus a liquidity premium. The liquidity premium can be viewed as an instantaneous price impact 
that affects transaction prices but not the value of the players' inventory. This form of market 
impact modeling is analytically more tractable than that of Obizhaeva & Wang (2006) which 
also allows for temporary price impacts and resilience effects. It has also been adopted by, e.g., 
Carlin, Lobo & Viswanathan (2007), Schoeneborn & Schied (2008) and many practitioners from 
the financial industry. 

Our framework is flexible enough to allow for rather general liquidity costs including the linear 
cost function of Almgren & Chriss (2001) and some form of spread crossing costs. We show 
that when the market participants are risk neutral or have CARA utility functions the pricing 
game has a unique Nash equilibrium in the class of absolutely continuous trading strategies; 
existence results for more general utility functions are given for the one player case. We solve 
the problem of equilibrium pricing using techniques from the theory of stochastic optimal control 
and stochastic differential games. We show that the family of the players' value functions can 
be characterized as the solution to a coupled system of non-linear PDEs. Here we use a-priori 
estimates for Nash equilibria; we prove that the system of PDEs has a unique classical solution 
with bounded derivatives. It turns out that the equilibrium problem can be solved in closed form 
for a specific market environment, namely the linear cost structure used in Almgren & Chriss 

(2001) and risk neutral agents. We use this explicit solutions to state some conditions which make 
manipulation unattractive or avoid it altogether. For instance, we show that when the agents 
are risk neutral no market manipulation occurs in zero sum games, i.e., in a game between an 
option writer and an option issuer. Furthermore, we find that the bid ask spread is important 
determinant of market manipulation. It turns out that the higher the spread, the less beneficial 
market manipulation: high spread crossing costs make trading more costly and hence discourage 
frequent re-balancing of portfolio positions. 
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This paper is organized as follows: We present the market model in section 2. In chapter 3, 
we formulate the optimization problem, derive a priori estimates for Nash equilibria and prove the 
existence of a solution for one player with general utility function. We solve the multi-player case 
in section 4 for risk neutral and CARA agents. We use these solutions in chapter 5 to show how 
market manipulation can be reduced. Section 6 concludes. 

2 The Model 

We adopt the market impact model of Schoeneborn & Schied (2007) with a finite set J of agents, 
or players, trading a single stock whose price process depends on the agent's trading strategies. 
Following Almgren & Chriss (2001) we shall assume that the players have a permanent impact on 
asset prices and that all trades are settled at prevailing market prices plus a liquidity premium 
which depends on the change in the players' portfolios. In order to be able to capture changes 
in portfolio positions in an analytically tractable way, we follow Almgren & Chriss (2001) and 
Schoeneborn Sz Schied (2007) and restrict ourselves to absolutely continuous trading strategies. 
Hence we consider only trading strategies from the class 

X = {X : [0, T] h->- R\X absolutely continuous, adapted and Xq = 0} 

where we denote by Xf the number of stock shares held by player j £ J = {1,...,N} at time 
t G [0, T], We write dX\ = X\dt and call X^ the trading speed of the player j. 

2.1 Price dynamics and the liquidity premium 

Our focus is on valuation schemes for derivatives with short maturities under strategic market 
interactions. For short trading periods it is appropriate to model the fundamental stock price, i.e., 
the value of the stock in the absence of any market impact, as a Brownian Motion with volatility 
(aBt). Market impact is accounted for by assuming that the investors' accumulated stock holdings 
Y^iLi X 1 have a linear impact on the stock process (Pt) so that 

N 

Pt = Po + aB t + \Y J Xl (2-1) 

i=i 

with a permanent impact parameter A > 0. The linear permanent impact is consistent with the 
work of Huberman & Stanzl (2004) who argued that linearity of the permanent price impact is 
important to exclude quasi-arbitrag^J 

2 There is some empirical evidence that very large trades have a concave price impact but this observation needs 
further validation. 
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A trade at time t G [0, T) is settled at a transaction price Pt that includes an additional 
instantaneous price impact, or liquidity premium. Specifically, 



Pt=p t +g(y j xi\ (2.2) 



with a cost function g that depends on the instantaneous change Yli=\ X 1 in the agents' position 
in a possibly non-linear manner. The liquidity premium accounts for limited available liquidity, 
transaction costs, fees, spread crossing costs, etc. Spread-crossing costs are of particular impor- 
tance and have not been considered in the previous literature on market impact. 

Remark 2.1. In our model the liquidity costs are the same for all traders and depend only on the 
aggregate demand throughout the entire set of agents. This captures situations where the agents 
trade through a market maker or clearing house that reduces the trading costs by collecting all orders 
and matching incoming demand and supply prior to settling the outstanding balance Yli=i ^-t °^ 
market prices. o 

We assume with no loss of generality that g is normalized, i.e., g(0) = and smooth. The 
following additional mild assumptions on g will guarantee that the equilibrium pricing problem 
has a solution for risk neutral and CAR A investors. 

Assumption 2.2. • The derivative g' is bounded away from zero, that is g' > e > 0. 
• The mapping z \— > g(z) + zg'(z) is strictly increasing. 

The first condition is natural for a cost function. Since the liquidity costs associated with a net 
change in the overall position z is given by zg(z), the second assumption states that the agents 
face increasing marginal costs of trading. Our assumptions on g are satisfied for the following 
important examples: 



Example 2.3. Cost functions which satisfy Assumption 2.2 are the linear cost function g{z) = nz 



with k > 0, used in Almgren & Chriss (2001) and cost functions of the form 

2 

g{z) = kz + s— arctan(Cz) with s, C > 0. 

7T 

The former is the cost function associated with a block-shaped limit order book. The latter can 
be viewed as a smooth approximation of the map z i— > kz + s ■ sign(z) which is the cost function 
associated with a block-shaped limit order book and spread s > 0. 

2.2 Preferences and endowments 

Each agent is initially endowed with a contingent claim W = H^Pt), whose payoff depends on 
the stock price Pt at maturity. Although it is not always necessary we assume that the functions 
H J are smooth and bounded with bounded derivatives Hp. 
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Remark 2.4. We only consider options with cash delivery. The assumption of cash delivery is 
key. While cash settlement is susceptible to market manipulation, we show in Section 5 below 
that when deals are settled physically, i.e., when the option issuer delivers the underlying, market 
manipulation is not beneficial: the cost of acquiring at an increased price outweighs the benefits 
from a possible higher option payoff, due to an increase in the underlying. o 

We model the risk preference of the agent j £ J with a von Neumann - Morgenstern utility 
function u 3 and assume that her aim is to maximize her expected portfolio value at maturity from 
trading in the financial market so the agent's optimization problem is given by: 

cT 



sup IE 



X\P t dt + H 3 (P T ) + value 3 (X 3 ,) 



(2.3) 



The portfolio value consists of the trading costs — J Q X 3 P t dt, the option payoff H 3 (Pt) and 
the liquidation value value 1 (X° T ) of X 3 T stock shares at maturity. Rigorously defining a form 
of liquidation value in a multi-player framework is a question of its own mathematical interest, 
and is not the focus of this paper. The problem of defining a liquidation value in a single-player 
framework has been solved in a recent paper by Schoeneborn & Schied (2008); in a game-theoretic 
setting the problem is much more involved. The agents optimize against their beliefs about the 
other players' individual assessments of their respective portfolio values. In order to simplify the 



analysis we work under Assumption 2.5 below. 



Assumption 2.5. All agents optimize their utility assuming that 



value 3 {XL 



T 



X 3 P t dt. 



(2.4) 



The preceding assumption is motivated by the single player, risk-neutral framework where (2.4) 



holds in expected values. In such a setting it states that the expected costs (utility) of building 
up a portfolio over the time span [0,T] under market impact equals the expected liquidation 



costs (utility) under infinitely slow liquidation. In a multi-player model Assumption 2.5 is only 
a first benchmark that simplifies the subsequent analysis even if all player are risk-neutral. The 
assumption is nonetheless useful as it allows us to carry out our analysis of equilibrium and to 
derive some insight into the structure of market dynamics under strategic interactions. 



Under Assumption 2.5 the individual optimization problem (2.3) reduces to 



sup E 

X3£X 



£x 3 g dt + H 3 (P T ) 



(2.5) 



We say that a vector of strategies (X 1 , ...,X N ) is a Nash equilibrium if for each agent j G J 
her trading strategy X 3 is a best response against the behavior of all the other players, i.e., if 
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Xi solves (2.5), given the vector X~i = (X 1 )^. In the following section we derive a-priori 
estimates for equilibrium trading strategies and use standard results from the theory of stochastic 
optimization to show that Nash equilibria can be characterized in terms of a coupled system of 
partial differential equations (PDEs). For the special case of risk neutral and CARA investors 
we show that the system of PDEs has a solution so that a unique (in a certain class) equilibrium 
exists. 



3 Equilibrium Dynamics and A-Priori Estimates 



In this section we formulate the optimization problem (2.5) as a stochastic control problem, derive 
the associated Hamilton-Jacobi-Bellman-equations, HJB for short, and transform it into a system 
of coupled PDEs. To this end, we choose the stock price P and the trading costs W of the agent 
j £ J as state variables. They evolve according to: 

N 

dP t = adBt + Xj^Xidt, P =Po (3.1) 



»=i 



(3.2) 



For a given time t < T, spot price p and trading costs r the value function of the player j, defined 

by 



V J (t,p, r) = sup K t 



[ X s9 (E *lj ds + &(P T )\ \P t =p 



(3.3) 



denotes the maximal expected portfolio value at maturity that the player can achieve by trading 
the underlying. The associated HJB-equation is 

1 



= v t + ^ 2 «pp + sup 



A (c j + X~ j vi + (?q (9 + X~ j v J 



v J (T,p,r) = v? (— r 3 +H J (p)^ 



(3.4) 



cf. Fleming & Soner (1993). Here we have used v 3 i = for % ^ j. Given the trading strategies 
X~i of all the other agents, a candidate for the maximizer e 7 = X 3 should satisfy 



= X-Z- + g ( d> + X- j ) + cV ( c? : + X- 



(3.5) 



provided v J rJ / 0. We sum up these equations over the set of players in order to get the following 



characterization for the cumulated equilibrium trading speed J2iLi X't- 



N i 
V 



N 



N 



N 



i=l r 



(3.6) 



vi=l 
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4=1 ~r \ • Plugging this solution back into (|3.5|) allows to compute 

r 

the optimal strategy c 3 = X 3 in terms of v l p and as 



Due to Assumption 2.2 z ^ Ng{z) + zg'(z) is strictly increasing. Hence, equa tion (3.6) admits a 
unique solution X* 



c 3 =X 3 



1 



g'lx 



(3.7) 



This expression is well defined if v a / because g' > 0. To conclude, we have turned the family 



of individual HJB-equations (3.4) into the following system of coupled PDEs for j = 1, N: 

1 







vi + -o~ 2 vi 



} !rp + \X*v p + X 3 g[X> 
^ v 3 (T,p, r) = u 3 (—r 3 + H 3 (p)) 



where the coupling stems from the aggregate trading speed X* = X* ( ^ 



(3.8) 



Solving the system (3.8) is delicate, to say the least. The problem is the non-linearity coming 
from the expression X 3 g along with the implicit dependence of X 3 on the derivatives v l p 

and v l i and the fact that v l ri appears in the denominator. The latter problem can be coped with 
by assuming that (after a possible monotone transformation) the agents' risk preferences satisfy 
a translation property so that u r = 1. A large class of such utility functions can be linked to 
backward stochastic differential equations but we choose not to embed our work into that line of 
research. We consider instead the case of risk neutral and CARA investors where the existence 



of a unique classical solution to the system (3.8 ) is guaranteed without any reference to backward 



equations. The proof uses the following a-priori estimates for the optimal trading strategies. It 
states that in equilibrium the trading speed is bounded. As a result, the agents' utilities from 
trading and the value function associated with their respective HJB equations along with (as we 
shall see) their derivatives are bounded. 



Lemma 3.1. 

satisfies 



Let (^X 1 , X N ^j be a Nash equilibrium for problem (2.5). Then each strategy X 3 



X 3 



< N— max \ \H, 



£ 



Plloo ' 



Proof. Let j £ J, h = maxj 1 1 iJ * [ | ^ and A = j^C^i^* — o} t> e the set where the aggregate 
trading speed is nonnegative. Let us fix the sum of the competitors' strategies X~ 3 . On the 
set A the best response X 3 is bounded above by K = jh. Otherwise the truncated strategy 
Y 3 = X 3 A jh would outperform X 3 . To see this, let us compare the payoffs associated with X 3 
and Y 3 . The payoff associated with X minus the payoff associated with Y can estimated from 
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below as 



Y j g [Y j + X~ j ) dt + H j (P T (Y j )) 



o 

+ J X j g [x j + X- j ^j dt - W(P T (X j )) 

> / Y j (g [X j + X~A - g (t j + X~ j 

+ (X* - Y 3 ) g (X* + X-*) dt - \(X 3 T - YJ>) . 

Note that Xi + X~i > on A and thus g (x* + X~^j > due to Assumption 

g (X* + X~A - g (iri + X~A > e fx* - Y^ , again by Assumption 

payoffs is therefore strictly bigger than 

rT 



2.2 



2.2 



Furthermore, 



The difference in the 



o 



Y J e (X 3 -Y 3 \dt 



Xh I X J - YHt 
'o 



> 



Xi>Yi 



:Y J - A/i) (X - Y')<II 



0. 



This shows that X 7 is bounded above on the set A. A symmetric argument shows that on the 
complement A c the optimal response is bounded below by —K. Furthermore, this implies that on 
A c the optimal response is bounded above by 

N 

x j = j2 xi + - + ( N ~ 

i=l i^j 

A similar argument yields that X 7 is bounded below on A. This completes the proof. □ 
Remark 3.2. These a priori estimates, together with the boundedness assumption on the payoffs 



H 3 , imply that the optimization problem (2.5) is bounded. Moreover, if a smooth solution to the 
PDE system (3.8) exists, it is bounded. o 



In the one player framework we can use a standard result from the theory of stochastic control 

= u{-r + H(p)) letip e <T 3 



to show that (3.8) admits a unique solution. 

1 and for the terminal condition tp(p, r 



Proposition 3.3. Let N 

and let ip, ip p , tp r satisfy a polynomial growth condition. Then the Cauchy problem \3.8^ admits a 
unique classical solution in "Jf ' 2 , which coincides with the value function V. 



Proof. Due to Lemma (3.1), the optimal control X can be chosen from a compact set. Thus, we 



can apply Theorem JV.4.3 in Fleming & Soner (1993), which yields that (3.8) admits a unique 
solution in If > 2 , which is of polynomial growth. It remains to apply the Verification Theorem 
IV.3.1 from Fleming & Soner (1993) to see that this solution coincides with the agent's value 
function V. □ 
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4 Examples 



In this section we establish existence and uniqueness of equilibrium results for risk neutral and 
CARA investors. For risk neutral investors and linear cost functions the equilibrium strategies can 
be given in closed form; if spread crossing costs are involved a closed form solution is not available 
and we report numerical results instead. 



4.1 Risk Neutral Agents 

Let us assume that all players are risk neutral, i.e. u J (z) = z. In a first step we prove existence of 



a unique solution to the system (3.8) for general cost functions g. Subsequently we construct an 



explicit solution to (3.8 1 for the linear cost structure used in Algmen & Chriss (2001). 



4.1.1 General Cost Structure 



Let g be a general cost function which satisfies Assumption |2.2| In the risk neutral case the value 
function of player j turns into 



V J (t,p, r) = — r 3 + sup Et 



[ ±j s9 (E *lj ds + H\P T )\P t = p 



In particular, = —1 so the optimal strategies do not depend on the trading costs. In other 
words, the state variable r is redundant and we omit it in this section. We write 



The HJB-equation (3.4) turns into 

= 4 + \ a2v l P + SU P 



A ( cP + X-3 ) v 3 p - c 3 g ( c? + X-i 



where we have used v J rj = —1. The optimal trading speed form (|3.7|) is given by 



d> = X j 



1 



g'(X* 

N 



-H + g(X*) 



where the aggregate trading speed X* = ^2 i=1 X 3 is the unique solution to 



jv / N \/ N \/ N 

i=l \i=l ) \i=l I \i=l 



XI 



The system of PDEs (3.8) therefore takes the form 

1 



= vj + -a 2 vi„ + XX* vi, - X j g X* 



J t 1 2" "pp 1 "v 



(4.1) 



(4.2) 



(4.3) 



(4.4) 
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with terminal condition v°{T,p) = W{p). The following proposition shows that a unique solution 
exists if H 6 i.e. H and its derivatives up to order 2 are bounded. The proof follows from a 
general existence result stated in Appendix A. 



Proposition 4.1. Let H E "jf 6 2 . T/ien the Cauchy problem (4-4) admits a unique classical solution 
in c € 1 ' 2 , which coincides with the vector of value functions. 



An alternative way of solving the system (4.4 1 is the following: If we sum up the N equations, 



we get a Cauchy problem for the aggregate value function v = Y^i=i v% i namely 

= v t + -a 2 Vpp + X* [\v p - g (1*) 

with terminal condition v(T,p) = J2iLi^(p)- Existence and uniqueness of a solution to this 
one-dimensional problem can be shown using Theorem IV. 8.1 in Ladyzenskaja (1968). Once the 



solution is known, we can plug it back into (4.4) and get N decoupled equations. This technique 



is applied in the following section where we construct an explicit solution for linear cost functions. 



4.1.2 Linear Cost Structure 

For the particular choice g(z) = kz (k > 0) used in Algmren & Chriss (2001) and Schoeneborn 



Schied (2007), the solution to (4.4) can be given explicitly. 



Corollary 4.2. Let g(z) = kz. Then the solution of (4-4) can be given in closed form as the 
solution to a nonhomogeneous heat equation. 



Proof. The optimal trading speed from \A.2\ and the aggregate trading speed from (|4.3|) are 

X j = 



X* 




N 



(4.5) 



A 



N 



i=l 



k{N + 



A 



i=i 



k{N + 1) 



(4.6) 



Let us sum up the N equations from \AA\. This yields the following PDE for the aggregate value 
function v 



H=X 



v : 







1 2 X 2 N 2 



(4.7) 



with terminal condition v(T,p) = Yli=i^ z (p)- This PDE is a variant of Burgers' equation, cf. 



Rosencrans (1972). It allows for an explicit solution, which we cite in Lemma 4.3 With this 
solution at hand, we can solve for each single investor's value function. We plug the solution v 



back into the equations (4.5) and (4.6) for the trading speeds, and those into the PDE (4.4). This 
yields 







1 



ii 



with terminal condition v 3 (T,p) = H 3 {p). This nonhomogeneous heat equation is solved by 

A 2 " 



v 3 (T-t,p) = / H J dJ^(p,at) + 



Vp(s, -)d^Y(p, a(t — s)) 



where v is given in Lemma 4.3 and JV denotes the heat kernel. 



□ 



In the preceding proof and in Corollary 4Jj we need the solution to a variant of Burgers' 
equation. We cite it in the following Lemma. 

Lemma 4.3. Let A > 0,B / 0. The PDE 



with terminal value 



is solved by 



v(t,p) 



A 
B 



log 



= 2v t + Av pp + Bv, 



v(T,p) = G(p) 



,T-t 



Proof. By means of a linear transformation we can reduce the problem to A = B = 1. This 
particular case is solved in Rosencrans (1972). □ 



4.1.3 Numerical Illustrations 

In the risk neutral setting, we were able to reduce the system of PDEs from the multi-player setting 



to the one-dimensional PDE (4.7) for the aggregate value function. This can be interpreted as the 



value function of the representative agent. Such reduction to a representative agent is not always 
possible for more general utility functions. In the sequel we illustrate the optimal trading speed 
X(t,p) and surplus^] of a representative agent as functions of time and spot prices for a European 
call options H{Pt) = (Pt — K) + and digital option H(Pt) = ^{p t >k}, respectively. We choose 
a linear cost function, strike K = 100, maturity T = 1, volatility a = 1 and liquidity parameters 
A = k = 0.01. We see from Figure [T] that for the case of a call option both the optimal trading 
speed and the surplus increases with the spot; the latter also increases with the time to maturity. 
Furthermore, the increase in the trading speed if maximal, when the option is at the money. For 
digital options the trading speed is highest for at the money options close to maturity as the trader 
tries to push the sport above the strike. If the spot is far away from the strike, the trading speed 
is very small as it is unlikely that the trader can push the sport above the strike before expiry. 

3 By surplus we difference between the representative agent's optimal expected utility v(t,p) and the conditional 
expected payoff Kt[H(Pr)\Pt = p] in the absence of any market impact. 
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Trading Speed 



Surplus 




Figure 1: Trading speed and surplus for a risk neutral investor holding a European Call option. 



Trading Speed 




Figure 2: Trading speed and surplus for a risk neutral investor holding a Digital option. 
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Trading Speed at t=0 Surplus at t=0 




SPOT 



Figure 3: Trading speed and surplus for a risk neutral investor holding a European Call option for different 
spread sizes s = (black), 0.001 (blue), 0.002 (red), 0.003 (green), 0.004 (brown). The higher the spread, 
the smaller the trading speed and the surplus. 



Figures [3] and [4] illustrate that a high spread makes manipulation unattractive. It shows the 
optimal trading speed and the surplus at time t = for the Call and Digital option in the one 
player framework. We used the cost function 

g(z) = kz + s- sign(z) for different spreads s G {0, 0.001, 0.002, 0.003, 0.004} 

with the remaining parameters as above. We see that the higher the spread, the smaller the 
trading speed and the surplus. This is intuitive as frequent trading, in particular, when the option 
is at the money, incurs high spread crossing costs. The same is true for fixed transaction costs 
which also discourage frequent trading. 

4.2 Risk Averse Agents 

A second class which yields explicit results is those of exponential utility functions v?(z) = 
— exp (— a? z) for j = 1,...,N, where a- 7 > is the risk aversion coefficient. In this case the 
value functions satisfy 

V j {t,p,r) = exp (aV) • V j (t,p,0) 

and thus = a^V^. We suppress the state variable r and write V^{t,p) = V^(t,p, 0). As above, 
we first show existence and uniqueness of a solution for a general cost structure. In a second step, 
we derive the closed form solution for the linear cost function in the single player framework. 
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Trading Speed at t=0 Surplus at t=0 




SPOT 



Figure 4: Trading speed and surplus for a risk neutral investor holding a Digital option for different spread 
sizes s = (black), 0.001 (blue), 0.002 (red), 0.003 (green), 0.004 (brown). The higher the spread, the 
smaller the trading speed and the surplus. 



4.2.1 General Cost Structure 



The HJB-equation (3.4) turns into 

1 



= v 3 t + ^ 2 v j pp + sup A I c? + X~ j J v J p + cP g ( c? + X' j J a j v 3 



(4.8) 



with terminal condition v^(T,p) = — exp (— a J '.fP(p)). We apply the linear transformation V 
1 log(— ip) to turn the HJB equation into 



a? 



1 I 

u 

r pp 2 



= v> t + ^a 2 vi„ - ^crV (vV 



c3 



) + sup A [cP + X~ j ) v> p - c J g [d> + X- j 



(4.9) 



with terminal condition v 3 (T,p) = H 3 (p). Note that this equation equals the H JB- equation (4.1) 

2 



in the risk neutral setting, up to the quadratic term —^a 2 a 3 yJpj ■ As in (4.2), the optimal 
trading speeds are 



d> = X j 



1 



-Xvi + g(X*) 



9' [X^ 

where the aggregate trading speed X* is the unique solution to 



(4.10) 



N 



N 



N 



N 



o = xJ2v;-n 9 - K E^ 



(4.11) 



1=1 



vt=l 



vi=l 



vi=l 



If we plug X* and X- 7 back into (4.9), we get 

1 o_, 1 



= v\ + ^cr 2 ^ - -aV + AX*u£ - Xi g (X* 



2 pp 2 

We can show existence and uniqueness of a solution. 



(4.12) 
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Proposition 4.4. Let ffl £ ^ /or eac/i j G J. The Cauchy problem 4-9 admits a unique solution, 
which coincides with the vector of value functions (up to an exponential transformation). 



Proof. See appendix A. □ 
4.2.2 Linear Cost Structure, Single Player 

For the one player case with linear cost structure, we have an explicit solution: 



Corollary 4.5. Let N = 1 and g(z) = kz. Then the Cauchy problem 4-9 admits a unique solution, 
which can be given in closed form. 



Proof. The maximizer in (4.9) is 



c = X = ±* p 



and the Cauchy problem (4.12) turns into 



1 9 ~ ( 1 ■ > \ _ > 



= *>t + -o- v pp + I -a a I v p 



with terminal condition v(T,p) = H(p). This is Burgers' equation. Its explicit solution is given 
in Lemma l4~3l □ 



4.2.3 Numerical Illustrations 



Let us conclude this section with numerical illustrations. We simulated the system (4.8) for two 
players. Figure [5] shows the aggregate optimal trading speed X(0,p) + Y(0,p) and the surpluses 
vi(p, p) — E [-u- 7 (H(Pt)) \Pq = p] for time t = and different spot prices p G [95,105] for the 
European Call option H(Pt) = (Pt — K) + ; we assume that Player 1 (blue) is the option writer 
and Player 2 (red) the option issuer. We chose the strike K = 100, maturity T = 1, volatility 
a = 2 and liquidity parameters A = k = 0.01 and risk aversion parameters a 1 = 0.01, a 2 = 0.01 
(plain), respectively, a 1 = 0.001, a 2 = 0.1 (dashed). Since Player 1 has a long position in the 
option, she has an incentive to buy the underlying; for the same reason Players 2 has an incentive 
to sell it (Panel (b)). Our simulations suggest that the dependence of the equilibrium trading 
speed on the agents' risk aversion is weak (Panels (b) and (c)) and that overall the option issuer 
is slightly more active than the option writer. Furthermore, we see from Panel (d) that the issuer 
benefits more from reducing her loss than the writer benefits from increasing her gains. This effect 
is due to the concavity of the utility function and increases with the risk aversion. 
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Value Function at t=0 



Trading Speed at t=0 




Figure 5: Value function, trading speed, aggregate trading speed and surplus for the writer (blue) and 
issuer (red) of a European Call option when both agents are risk averse. The plain (dashed) curves display 
the case where issuer is about as (more) risk averse than the option writer. 
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5 How to Avoid Manipulation 



In this section, we use the closed form solutions for risk neutral agents derived in subsection |4 . 1 . 2 
to illustrate how an option issuer may prevent other market participants from trading against her 
by using their impact on the dynamics of the underlying. Some of our observations were already 
made in Gallmeyer & Seppi (2000) in a three-period binomial model. We start with the simplest 
case of a zero-sum game. 

Corollary 5.1. Let all players be risk neutral with offsetting payoffs Y2d = iH l = 0. Then the 
aggregate trading speed is Yli=i X 1 = 0. 



Proof. If J2iLi H l = then (|4.6[) and (|4.7[) imply that the aggregate value function v = J2 



N 



equals zero. It follows from (4.6) that the aggregate trading speed vanishes. □ 



Loosely speaking in a zero-sum game, if all option traders are risk neutral and willing to move 
the market in their favor, their combined effect cancels. We note that this is no longer true for 
general utility functions, as illustrated in figure [5j Of course, in reality market manipulation is 
illegal and many investors are unable to manipulate the underlying in the first place. This is why 
we now look at the following asymmetric situation: The option issuer, Player 0, does not trade the 
underlying; her competitor, Player 1, owns the derivative H 1 / and intends to move the stock 
price to her favor. In addition, there are N — 1 informed investors without option endowment 
in the market. They are "predators" that may supply liquidity and thus reduce the first player's 
market impact, cf. Carlin, Lobo Sz Viswanathan (2007) and Schoeneborn Sz Schied (2007). We 
find that the more informed competitors are active, the less aggregate manipulation will occur. 

Corollary 5.2. Let H 1 / and H* = for i = 2, N. Then lim^oo E*Li M = 0. 



Proof. The solution v to the Cauchy problem (4.7) with terminal condition 

N 



,{T,p) = Y,H i {p) = H\p) 



i=l 



is bounded uniformly in N. Equation (4.6) yields 

N 

kN + 1 



E<ri ^ 1 JV->oo _ 

X = — — — -v v -> 0. 
k TV 4- 1 p 



i=l 

□ 

Let us modify the preceding setting a little. Again, Player issues a product H and does not 
intend to manipulate the underlying, while her competitors do. More precisely, assume that player 
splits the product H into pieces and sells them to N risk neutral competitors, such that each of 
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cumulated Trading Speed at t=0 cumulated Trading Speed at t=0 




SPOT 



(a) Call (b) Digital 



Figure 6: Aggregate trading speed X* at time t = for N = 1 (black), 10 (blue), 100 (red) players each 
holding 1/N shares of a Call (left) and Digital (right) option with strike K = 100. The more agents, the 
less aggregate manipulation. 



them gets jjH. We find that their aggregate trading speed YliLi X % is decreasing in the number 
of competitors N. Consequently, the option issuer should sell her product to as many investors as 
possible in order to avoid being outsmarted. We illustrate this result in figure [6] which shows the 
aggregate trading speed at time t = of N players each holding 1/N option shares. 

Corollary 5.3. Let H % = for i = 1,...,N, Then for each t G (0, T] the aggregate trading 
speed YliLi is decreasing in N and limjv^oo Yli=i -^-t = 0- 



Proof. As in the preceding corollary, the assertion follows from equations (4.7) and (4.6). □ 



The preceding results indicate how an option issuer can prevent her competitors from manip- 
ulation. One strategy is public announcement of the transaction: the more informed liquidity 
suppliers on the market, the smaller the impact on the underlying. A second strategy is split- 
ting the product into pieces; the more option writers, the less manipulation. Let us conclude this 
section with a surprisingly simple way to avoid manipulation: using options with physical delivery. 

Remark 5.4. Calls, Puts and Forwards with physical delivery do not induce stock price manipu- 
lation, o 

Proof. Consider a risk neutral agent who owns > Call options with physical delivery and 
strike K. As above, we denote by X her strategy in the underlying. At maturity, she exercises 
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< 9 < 6 of her Call options. Problem (2.5) turns into: 

f 

/o 



sup E 

x&x,e<e 



X t {P t + nX t )dt + (X T + 0)[P t - -X(X t + 



9K 



sup E 

xex,e<e 



nXUt + 6[P T 



o 



1 



A0 





r f T 






sup E 




+ sup 


['( 


xex 


io 


6»<e 





-xe 



9K 
- 9K 



sup 



e\p T - -xe 



OK 



where the first term in the first line describes the expected trading costs in [0, T] and the liquidation 
value of 9 + Xt stock shares at maturity. The optimal trading strategy is X = 0. The same holds 
true for Put options and Forward options with physical delivery. □ 

6 Conclusion 

We investigated the strategic behavior of option holders in illiquid markets. If trading the under- 
lying has a permanent impact on the stock price, the possession of derivatives with cash delivery 
may induce market manipulation. We showed the existence and uniqueness of optimal trading 
strategies in continuous time and for a general cost function; in the one player framework for 
general utility functions, and in the multi-player case for risk neutral as well as CARA investors. 
Moreover, we showed how market manipulation can be reduced. 

Our work may be extended in several directions. Foremost, we derived our results under 



Assumption |2.5[ This assumption is only satisfied in the single-player risk-neutral case where the 
expected costs of buying a portfolio over a finite time interval under market impact equals its 
expected liquidation value under infinitely slow liquidation and does not hold in general. The 
problem of defining a proper notion of liquidation value under strategic interaction is important 
but was not our focus and is left for future research. Furthermore, it would be interesting to 
consider American or path-dependent options as well as more sophisticated market impact models 
such as Obizhaeva Sz Wang (2006) that account for resilience effects and temporary price impacts. 



A An Existence Result 



In this section, we prove Propositions 4.1 and 4.4 where the PDE (4.4) in the risk neutral setting 
is a special case of the system (4.9) for risk averse agents, with a- 7 = for each j. In order to 
establish our existence and uniqueness of equilibrium result, we adopt the proof of Proposition 
15.1.1 in Taylor (1997) to our framework. After time inversion from t to T — t both systems of 
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PDEs are of the form 



v t = Lv + F(v p ) 



(A.l) 



for v 



^ ( 1 



), where L is the Laplace-operator 



and F = (F 1 , F ) is of the form 



|aV(^) 2 + AXX-^ (**)• 



Here X* and X? are given implicitly by (4.2) and (4.3). The initial condition is 

v(0, P ) = H(p) = (H\...,H N ). 



We rewrite (A.l) in terms of an integral equation as 

^ L F(v p (s))ds = Vv(t). 



v(t) 



e tL + 



(A.2) 



(A.3) 



and seek a fixed point of the operator ^ on the following set of functions: 
X = rf(R,R N ) = {ve ^\R,R N ) I v,v p bounded} 
equipped with the norm 

IMIx - IMloo + Klloo • 

We set Y = Note that X and Y are Banach spaces and the semi-group e tL associated with the 
Laplace operator is strongly continuous on X, sends Y on X and satisfies 

II NjSf(Y,X) - 

for some C > 0, 7 < 1 and t < 1. Furthermore, the nonlinearity F is locally Lipschitz and 
belongs to < ^ 00 . Indeed, the map a \-> X*(a) is due to the implicit function theorem with 
first derivative 



_d_ 

dv r 



X*(v p ) 



(N + l)g'(X*(v p ) - X*(v p )g"(X*(v p )) 



where the denominator is positive due to Assumption |2.2| The cost function g is by assump- 
tion. In particular, the assumptions of Proposition 15.1.1 in Taylor (1997) are satisfied. 



The a-priori estimates of Proposition 3.1 yield that, if a solution v to (A.l ) exists, it is bounded 



in the sense ||«|| x < K. Therefore, we define 

X K = {v G X I ||u|| x < K) 
and choosing K large enough we may assume that the initial condition satisfies H G Xk- 



We are now ready to prove existence and uniqueness of a solution to (A.3). Proposition 15.1.1 



in Taylor (1997) gives a solution for a small time horizon [0,r], with r > specified below. We 
apply his argument recursively to extend the solution to [0,T]. 
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Proposition A.l. There is r > such that for each n G No, the PDE (A. 3) with initial condition 



(A. 2) admits a unique classical, bounded solution in on the time horizon [0,nr AT]. This 



solution coincides with the value function. 

Proof. 1. For n = 0, there is nothing to prove. Pick n G N such that nr < T. By induction, 
we can assume that there is a solution G Hk on the time horizon [0,nr]. In particular, 
the initial condition for the next recursion step = t/ n )(nr) is in X^. 

2. Fix 5 > 0. We construct a short time solution on the following set of functions: 

z (n+i) a | w G ([ nT> ( n +l)r],X) | u(nr) = «(t) - /i (n) || x < 5 V*} . 

To this end, we first show that * : z( n+1 ) — ► z( n+1 ) is a contraction, if r > is chosen small 
enough. 

For this, let t\ be small enough such that for t < t\ and any v G Xj<- we have 



||e^-i x <^. 

Here we used that e tL is a continuous semigroup and \\v || x < .fT. In particular, for v = h^: 



1 



e tL h (n) - /i (n) 



1 r 
< -8. 

x ~ 2 



For v G Z^ n+1 \ the derivative v p is uniformly bounded in the sense ||t>plloo — 11^ llx + ^ — 
K + 5. Hence, we only evaluate F on compact sets. By assumption, F is locally Lipschitz. 
In particular, F is Lipschitz on compact sets. In other words, there is a constant K\ such 
that for any v, w G Z^ n+l "> we have 

\\F(v p ) - F(w p )\\ Y < Kt \\v - w\\ x 

This implies, for w = 

\\F(v p )\\ Y < F(t&>) +Jfi v-h™ 



± K 2 . 



This, together with the boundedness assumption on e* L , yields 



J (H,)I %(»))d!/ 



< t||e tL || sup \\F(v p (y))\l 

nr<y<t 

< t x ~^CK 2 . 
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8 

2CK 2 



1 

1-7 



This quantity is < \5 if t < t 2 = ( 
Finally, it follows that for v G we have 



- /i (n) 



< 



e' L /i (n) - h {n) 
5. 



+ 



,(*-»)£ 



F(v p (y))dy 



~ 2 2 



This shows that \1/ maps Z( n+1 ) into itself. 

It remains to show that \& is a contraction. Let v,w E Z( ra+1 ). Then 

rt 



[F(v p (y))-F(w p (y))}dy 



TIT 

tL\ 



< t\\e^\\ sup \\F{v p {y)) - F{w p {y))\\ Y 

riT<y<t 

< t 1 -" / CK 2 sup \\v(y) - w( 

m<y<t 



The quantity t x -^CK 2 is < \ if t < r 3 = 
in Z( n+1 \ if t is small in the sense 



1 

1-7 



This proofs that ^ is a contraction 



< r = min{Ti,T2,T 3 }. 
Note that the time step r does not depend on n. It is the same in every recursion step. 

3. It follows that ^ has a unique fix point v in Z^ n+l \ In other words, we constructed a 
function v 6 ^([nr,(n + l)r],X) = ^ x \nr,(n + l)r] which solves the PDE (0| with 
initial condition u(s) = /i ( - n - ) = i/ n )(nr) on the time interval [nr, (re + l)r]. 

This solution is actually in ^ 1,2 ((rrr, (n + l)r] x E, M^) , due to Proposition 15.1.2 in Taylor 
(1997). Furthermore, v is bounded by construction. Indeed, \\v {{^ < ||^^|| x + 5 < K + 5. 
We define the new solution as 



v (n+l) A v (n) lm 



t<riT} 



+ vt 



{riT<t<(n+l)T}- 



By construction, i;( n+1 ) solves (A.3) on the time horizon [0, (n + l)r] and is bounded and in 
< ^ 1,2 . Hence, we can apply the Verification Theorem IV. 3.1 from Fleming & Soner (1993), 
which yields that v^ n+1 ' coincides with the vector of value functions (up to time reversal and 
an exponential transformation, if a J > 0). The a priori estimate in Proposition 3.1 yields 
that r/ n+1 ) G X K . In particular, ||t>( n+1 )((n + l)r)|| x < K, which is necessary for the next 
recursion step. 

This completes the proof. □ 
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